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Abstract 

In flight control, the design objective and the aircraft dynamics may be different in low and 
high angle of attack regions. This paper presents a systematic switching LPV control design 
method to determine if it is practical to use for flight control designs over a wide angle of attack 
region. The approach is based on multiple parameter-dependent Lyapunov functions. A family 
of LPV controllers are designed, and each of them is suitable for a specific parameter subspace. 

The state of the controller is reset to guarantee the stability requirement of the Lyapunov 
function when the switching event occurs. Two parameter-dependent switching logics, hysteresis 
switching and switching with average dwell time, are examined. The proposed switching LPV 
control scheme is applied to an F-16 aircraft model with different design objectives and aircraft 
dynamics in low and high angle of attack regions. The nonlinear simulation results using both 
switching logics are compared. 
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1 Introduction 


In flight control, different performance goals are often desirable for different angle of attack regions. 
For example, in a low angle of attack scenario, pilots desire fast and accurate responses for maneu- 
vering and attitude tracking. While in a high angle of attack region, the flight control emphasis 
lies in the maintainability of aircraft stability with acceptable flying qualities. A modern fighter 
aircraft usually works in a wide angle of attack region, even near stall or poststall regime. In such 
a case, it is difficult to design a single robust controller over the entire flight envelope. Typically, 
the controller is designed by compromising the performance in some angle of attack region. 

Another issue encountered in flight control is that the actuator dynamics may be different in low 
and high angle of attack regions. It is known that the redundant control effectors, such as thrust 
vectoring nozzles, are usually incorporated in the high angle of attack region to provide additional 
control power. The usual way to generate the thrust vectoring command is a two-step procedure. 
A controller is designed first based on the generalized control, and the real control input is then 
generated using a control selector [1, 2, 3, 4], It has not been clearly addressed how to develop the 
control law with guaranteed stability and performance by considering the aerodynamic force and 
thrust force in a unified frame. 

To solve those problems, this paper considers a switching linear parameter-varying (LPV) control 
technique to design a family of controllers, each suitable in different angle of attack region, and 
switch among them according to the evolution of angle of attack. The whole framework is based 
on LPY systems because of their relevance to nonlinear systems. The proposed switching LPV 
control technique is the generalization of results in switched LTI systems. Obviously, stability 
is an important and challenging problem in switched systems, and it has received considerable 
attention in the recent literature (see Refs. [5, 6, 7]). For a family of stable LTI systems, the 
existence of a common Lyapunov function provides sufficient conditions for stability of switched 
systems under arbitrary switching sequences [8, 6]. However, this kind of stability guarantee is 
deemed to be too conservative when a particular switching logic is concerned. For restricted class 
of switching signals, multiple Lyapunov functions have been shown to be very useful tools for 
stability analysis purpose. Specifically, non-traditional stability conditions have been developed 
using either piecewise continuous Lyapunov functions [9, 10, 11, 12], or discontinuous Lyapunov 
functions [13, 14], In the former case, the values of Lyapunov functions corresponding to the active 
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subsystems form a decreasing sequence. This constraint is relaxed in the latter one by requiring 
each Lyapunov function V p to decrease when the pth subsystem is active with possible jump at 
switching time. The results of switched LTI systems have been generalized to the analysis and 
control of switched LPV systems [15], which is further extended in Ref. [16] by incorporating 
average dwell time switching logic [17]. The stability of switched LPY systems was analyzed by 
multiple parameter-dependent Lyapunov functions, which are allowed to be discontinuous at the 
switching surfaces. 

Recently, the authors have developed switching LPV control techniques under hysteresis and av- 
erage dwell time switching logics [18]. The multiple parameter-dependent Lyapunov functions [16] 
are used for analyzing the stability of switched LPV systems. We also considered the optimized 
performance, which is another important issue and has not been addressed adequately in the most 
of current work. However, the control synthesis conditions for both switching logics are formu- 
lated as non-convex matrix optimization problems, which are difficult to solve. In this research, 
we will bypass the non-convex condition by resetting the state of the controller when the switching 
event occurs. If the plant state is measurable, we can assign the controller state equal to the plant 
state at switching instances to guarantee the stability requirement of the Lyapunov functions. It is 
noted that the controller state reset is not always a feasible option due to lack of the information 
of the plant state. In this case, the state of the plant may be estimated by adding an observer. 
To simplify the theoretical derivation, we assume that the state of the plant is measurable in this 
research. The effect of an observer and the constraints for the observer design in order to guarantee 
the exponential stability will be studied in our future research. 

The notation is standard. R stands for the set of real numbers and R + for the non-negative real 
numbers. R TOXn is the set of real m x n matrices. The transpose of a real matrix M is denoted 
by M T . Ker(M) is used to denote the orthogonal complement of M . We use S nxn to denote the 
real symmetric n x n matrices and S” xn to denote positive definite matrices. If M G S nxra , then 
M > 0 (M > 0) indicates that M is positive definite (positive semidefinite) and M < 0 (M < 0) 
denotes a negative definite (negative semidefinite) matrix. For x G R", its norm is defined as 
||rr|| := (x T x)? . The space of square integrable functions is denoted by £ 2 , that is, for any u G £ 2 , 
IMI2 := [Jo°° u T (t)u(t)dt] 2 is finite. 

The paper is organized as follows: the problem addressed in this paper is first briefly formulated 
in Section 2. In Subsection 2.1, we provide the synthesis condition for hysteresis switching LPV 
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control via controller state reset, and the result for switching with average dwell time is provided 
in Subsection 2.2. The proposed technique is applied to an F-16 longitudinal flight control system 
with thrust vector augmentation in Section 3. The nonlinear simulations of the closed-loop system 
under these two switching logics are presented and also compared. Finally, the paper concludes in 
Section 4. 

2 Switching LPV Control via Controller State Reset 

Consider a generalized open-loop LPV system as functions of the scheduling parameter p. It 
is assumed that p is in a compact set V C R s with its parameter variation rate bounded by 
E.k < Pk < Pa- k = 1,2,--- ,s. The time-varying parameter value is assumed measurable in 
real-time. For notational purposes, we denote V = {v : < fk < Pjt, k = 1,2,..., s}, where V 

is a given convex polytope in R s that contains the origin. Suppose that the parameter set V 
is partitioned into a finite number of closed subsets {Vi} ieZN by means of a family of switching 
surfaces S tj (i,j E Z n), where the index set Zft = {1, 2, . . . , N}. In each parameter subset, the 
dynamic behavior of the system is governed by the equation 

x(t) A,j{p) B u {p) B 2 ,i(p) x(t) 

e{t) = C u (p) D lu {p) D 124 {p) d(t) , VpeVi , (1) 

y(t)_ C 2}i {p) D 2U (p) D 2U {p)_ «(*)_ 

where the plant state x E R n . e E R ne is the controlled output, and d E R n<i is the disturbance 
input, y E R"® is the measurement for control, and u E R n “ is the control input. All of the 
state-space data are continuous functions of the parameter p. Note that each LPV model should 
have the same number of state, and the reason will be clarified in the sequel. It is also assumed 
that 

(Al) (Aj(p), B 2 ^(p), C 2 ^(p)) triple is parameter-dependent stabilizable and detectable for all p. 

(A2) The matrix functions [B' 2 i ( p) Dj 2 y{p)\ and [C 2 y (p) D 2 \^(p)] have full row ranks for all 
p E Vi- 

(A3) Di U (p) = 0 and D 224 {p) = 0. 

Note that the assumption (A3) can be relaxed by loop transformation [19], but the controller 
formula will become more complicated. 
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Given the open- loop LPV system (1), it is sometimes hard to find single LPV controller working 
for the entire parameter region based on a single Lyapunov function (quadratic or parameter- 
dependent) [20, 21]. Switching LPV control technique permits use of most suitable controllers in 
different parameter subsets, and switch among them according to the evolution of the parame- 
ters. It will be beneficial to improve controlled performance and enhance design flexibility. In the 
previous research of switching LPV control [18], the state of the controller was preserved before 
and after switching, and this resulted in non-convex matrix optimization problems. In this paper, 
the controller state is reset at switching instances. As a result, this modification leads to a linear 
matrix inequality (LMI) optimization problem, which can be solved efficiently. The switched LPV 
controllers with state reset are in the form of 


Xk(t) 


AApA) 

Bk,i{p) 


Xk(t) 

u{t) 


Ck,i {p) 

V>k,i {p) 


y(t) 


[ Xk(t + ) = x(t), 


V p G Vi 


V p E S t j. 


( 2 ) 


where the dimension of controller state is Xk E R n * , and = n. Each controller is activated in a 
specific parameter subset V u and the state of the controller evolves continuously within the subset 
Vi- When the parameter trajectory hits the switching surface <Sij , the controller state is reset to 
the current value of the plant state. 


Then the closed-loop LPV system can be described by 



Xcl(t) 



. . 


{ x d {t + ) 

= [x{ 



Xcl(t ) 


d(t) 


\^cl,<j{PiP) B cUa (p) w _ 

, V P<E Vi 

CciAp) d ciAp)\ [ c A) J ( 3 ) 

); x(t)], V p E S; jr 

where x c i E R n + n * with x^ = [x T xj], and c is a switching signal, which is defined as a piecewise 

constant function. It is assumed that a is continuous from the right everywhere. The value 

of switching signal a represents the active parameter subset and thus determines the dynamic 
behavior of the plant and the controller. It is desirable that each controller stabilizes the open-loop 
system with best achievable performance in a specific parameter region, and meanwhile maintains 
the stability of the closed-loop system when switching the controller. 


A discontinuous Lyapunov function consisting of multiple parameter-dependent Lyapunov functions 
is used for stability analysis and control design of switched LPV systems. If there exist a family of 
positive definite matrix functions {Xi(p)} ieZn/1 and each of them is smooth over the corresponding 
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parameter subset V, . The multiple parameter-dependent Lyapunov functions can then be defined 
as 

V a (x,p) = x T X a (p)x. (4) 

Generally speaking, for a switched LPV system to be stable, the value of the discontinuous Lya- 
punov function V a is not necessarily decreasing over the entire parameter trajectory. In fact, it is 
often enough to require that the value of V a decreases in the active parameter region V, provided 
proper switching logic is adopted. This will lead to relaxed stability condition and provides en- 
hanced design flexibility. In this section, we consider the synthesis conditions of switching LPV 
control with two different switching logics by resetting the controller state. 

2.1 Hysteresis Switching 

When hysteresis switching logic is employed, it is assumed that any two adjacent parameter subsets 
are overlapped, as shown in Figure 1(a). Thus there are two switching surfaces between two adjacent 
parameter subsets. We use Sij to denote the switching surface specifying the one-directional move 
from subset V% to Vj. 



(a) Hysteresis switching region (b) Hysteresis switching signal 

Figure 1: Hysteresis switching region and switching signal a. 

The switching events occurs when the parameter trajectory hits one of the switching surfaces Sij or 
Sji . The evolution of the switching signal a is described as follows. Let <r(0) = i if p(0) E V , . For 
each t > 0, if a(t~) = i and p(t) E Vi, keep a(t) = i. On the other hand, if a(t~) = i but p(t) ^ Vi, 
i.e., hitting the switching surface Sij , let a(t) = j. Similarly, if a(t ) = j but p(t) ^ Vj, i.e., hitting 
the switching surface Sji, let a(t) = i. Obviously, the switching event is parameter-dependent, 
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as shown in Figure 1(b). For LPV systems, it is conceivable that parameter-dependent switching 
is more practical than state-dependent or time-dependent switching. Since a changes its value 
only after the continuous trajectory has passed through the intersection of adjacent subsets V, and 
Vj, chattering is avoided. Also due to bounded parameter variation rates, only finite number of 
switches will happen in any finite time interval. 

For the closed-loop system (3) under the hysteresis switching logic, if on the switching surface 
Sij , we have Vjfx c i,p) > Vj([x; x\, p), i.e. the Lyapunov function of the closed-loop system (3) 
non-increases when switching from V, to Vj. Then the jth controller can be activated safely. The 
synthesis condition of switched LPY control with hysteresis switching is given in the following 
theorem. 


Theorem 1 Given an open-loop LPV system. (1), the parameter setV and its overlapped partition 
{P i)iez N ’ if there exist positive definite matrix functions it!,;, 5,; : R s — > S" xn ,f € Zjy, such that for 
any p E Vi, 


K,(p) 


OR, 


Ri{p)Aj (p) + Ai{p)Ri{p ) - {Pk, "A Rfip^lfip) B u {p) 


k= 1 


C u {p)R,{p) 

bUp) 


~i I a „ o 

0 ~R^n c i 


MrApXO (5) 


A/1» 


dSi 


Aj (p)Si(p) + Si(p)Ai(p) + X, »k} Si(p)B u (p) Cjfip) 


k = 1 


bUp) Slip) 

CiAp) 


l''Jn d 

0 


0 


KsAp) < 0 (6) 


Hi In 

Ri{p) In 

In Si(p) 


>0, (7) 


where N R j(p) = Kei Bjj(p) D l%i(p) 0 and ^sAp) = Ker C 2 ,i{p) D 2 i A p) 0 , and for any 
P & <Sij , 

Rfip) < Rj (p) , (8) 

then the closed-loop LPV system. (3) is exponentially stabilized by switched LPV controllers with 
state reset in the entire parameter set V, and its performance is maintained as ||e||2 < 7 1 1 ci 1 1 2 with 
7 = max {7 i} ie ? n ■ 


Proof: The proof is similar to that in Ref. [18]. Therefore, we only show that the Lyapunov 
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function on the switching surface Sij satisfies Vj(x c i, p) > Vj([x\x], p) given the boundary condition 
(8). Following the idea in Ref. [18], we partition the Lyapunov function matrices of the closed- loop 
system (3) according to the dimensions of the plant and controller states as 


Xi(p) 


Slip) Ni{p) 
Nf(p) ? 


X~\p) 


Ri(p) M,j{p) 
mJ {p) ? 


where MAp)Nf{p) = I - RAp)Si(p), and “?” means the elements we don’t care. By choosing 
Mi{p) = Ri(p) and N, (p) = R~ 1 (p) — S', (p). the Lyapunov function Xi(p) can be written as 


Mp) = 

which can be further decomposed to 

Xi{p) = 


Si(p) Ri 1 {p)-S i {p) 
R-\p)-S l {p ) Si{p)-Rr\p) 


I -I 


1 

o 

l --H 



R i(p) 0 

o sap) r,; ] (p) 



I 0 


1 

1 — 1 


Therefore, 


P ) 


x T {x k - x) 1 


r-Hp) o 

0 S a (p)-R-Hp) 



X 


1 

1 

i 


( 9 ) 


Assume that the parameter trajectory hits the switching surface S^. Then the Lyapunov function 
before switching is 


VAx d ,p) = x r R i 1 (p)x + {x k - x) T ( SAp ) - R t l {p)) (x k - x). 


Combining with the coupling condition (7) gives 

VAxcup) > x' R, '(p)x. 


If the states of the controller are reset when switching, i.e., x k = x, then the Lyapunov function 
after switching is 

Vj([x;x],p) = x T Rj 1 (p)x. 


Therefore, Vj.(x c i,p ) > Vj([x; x\, p) is satisfied at the switching surface Sij if R~\p) > RJHp), 
which is equivalent to the boundary condition (8). Q.E.D. 
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Note that the hysteresis-based switching LPV control via controller state reset is formulated as 
a convex optimization problem, and can be solved using efficient LMI numerical algorithms [22]. 
It should be pointed out that the functional dependency of matrices R, and S, to parameters 
p k ,k = 1, . . . , s has to be chosen before solving the problem. The notation Ylk=i i—ki 7 A } 711 

(5)-(6) represents the combination of derivative terms in the form of when v k is taken as 

either u k or v k . Therefore each inequality means 2 s different LMIs which must be checked. 

After solving matrix functions R-ifp) and Sj(p), the gains of switched LPV controllers can be 
constructed using the following formula [23] 

Ak,i{p,p) = | A I(p) - - N i(p)^f- + S i(p) [Mp) 

+B^ l (p)F i (p) + Li{p)C 2 M] Mp) + -Si(p) [B u (p) + Li(p)D 2U (p)] 

7 i 


+~CUp) [C u (p) + D l2 M)Fi{p)] M~ T (p) (10) 

BMp) = N-\p)S i (p)L i {p ) (11) 

C k ,i{p) = Fj(p)Ri(p)M~ T (p) (12) 

D k ,i{p)= 0, (13) 

where the matrix functions F, (p) and Lj(p) defined as 


Fi(p) = - {d1 24 (p)Di2Ap)) 1 biBlMR-Hp) + D^(p)CiAp)} 

Li(p ) = - [ 7 iSrHp^dP) + b u(p) D 2u(p)] (D2u(p)F>lAP)y 1 • 

To comply with the proposed controller state reset scheme, we need to choose particular realizations 
of LPV controllers with Mj(p) = Ri{p) and N,(p) = R~ 1 {p) — Sj(p). 


Remark 1 For general specification of Mj(p) andNj(p), the synthesis conditions (5)-(8) still hold. 
However, the controller state on the switching surface is not simply reset to the plant state. From 
the structures of the Lyapunov function matrix and its inverse, one can get 


Mp) = 


Vie z N . 


Slip) Mp) 

| _N?(P) -Ny P )R t (p)M- T (p)_ 

Combined with the constraint Mi{p)Nf (p) = I — Ri(p)Si(p), Xjfp) can be decomposed and further 
simplified to 


Mp) 


I -Rf(p)Mi(p) 


1 

So 

“• i 

o 

i 


1 

o 

1 

1 

O 

1 


0 -N?(p)R t (p)M~ l (p)_ 


-M- T (p)R~ l (p) I 
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Therefore, on the switching surface Sij, the controller state needs to be reset to 

x k {t + ) = Mj (p)R~ 1 (p)x(t). 

2.2 Switching with Average Dwell Time 

If the overlapped region between two adjacent parameter subsets shrinks, it eventually becomes 
a single switching surface, as shown in Figure 2(a). Different from hysteresis switching, here Sij 
and Sji represent the same switching surface between subsets V, and Vj regardless which direction 
the parameter trajectory moves from. As shown in Figure 2(b), the switching signal changes 
its value right after the parameter hits the switching surface Sij no matter which direction the 
parameter evolves from. This usually requires the continuity of Lyapunov function across the 
switching surfaces. To relax the continuity requirement, we will consider the switching logic with 
average dwell time [17, 16], and only limited number of switchings is allowed between a finite time 
interval. 



(a) Switching region with dwell time (b) Switching signal with dwell time 

Figure 2: Switching region and switching signal a with dwell time. 

Denote N a (T , t) as the number of switchings among subsets V, on an interval (t, T). The switching 
signal a has average dwell time r a if there exist two positive numbers Nq and r a such that 

N a (T, t)<N 0 + VO <t<T, 

Ta 

where Nq is called the chatter bound. This idea relaxes the concept of dwell time, allowing the 
possibility of switching fast when necessary and then compensating for it by switching sufficiently 
slow later on. 
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Basically, all we need is the fact that there exists a positive constant p such that Vj(p) < pVi(p) 
when switching from V, to Vj [7]. Due to the interchangeability of i and j, the Lyapunov function 
on the switching surface must satisfy 

-Vj(p) < V;ip) < pVj (p) V p G Sij, (14) 

P 

where p > 1 . It allows the change of Lyapunov function by p times of its value before switching. 
As a consequence, the average switching frequency over a finite time interval is limited to A- to 
compensate for possible increase of Lyapunov functions. 

Different from hysteresis switching, the synthesis condition of switching LPV control with average 
dwell time cannot be convexified under some conservative assumptions if the control state is pre- 
served when switching [18]. The next theorem shows that the non-convex synthesis condition can 
be bypassed by resetting the control state. It is an important step for actual implementation, and 
we will give the detailed proof. 


Theorem 2 Given scalars Ao > 0, p, > 1, an open-loop LPV system, (1), the parameter set V and 
its partition {Vi} ieZN , if there exist positive definite matrix functions Ri,Si : R s — > S” x, Vi G Zjy, 
such that for any p G Vi 


KM 


KM 


Ri(p)Af(p)+Mp)Ri(p) 
-K M T 'k} + a 0 Ri(p) 

C u (p)Ri{p) 

KM 

Aj {p)Si(p) + Si(p)Ai(p) 

S e IS 1 - 

+ y~l {Pki v k} M- + Ao Si(p) 


RMCUp) B U { P ) 


T iAn e 

0 


k = 1 


dpk 
BTMSi(p) 

CiAp) 


7 An d 

0 


0 


MM < o (is) 


smbiAp) c(M 


MM < o (16) 


0 

— 7 An 

RM In 

In Si(p) 


> 0, (17) 


where MM = BjM D n M 0 and MAp) = Kei C 2 Ap) d 21 Ap) 0 , and for any 

P G S-ij 

-Rj{p) < RM < pRj(p), 

p 


li 


(18) 



then the closed-loop LPV system, (3) is exponentially stabilized by switched LPV controllers with 
state reset in the entire parameter set V for every switching signal a with average dwell time 


and its performance is maintained as ||e ||2 < 7 IMI 2 with 7 = max { 7 *} • 


Proof: We still assume that the Lyapunov function of the closed-loop system (3) is in the form of 
(4). If on the switching surface Sij, the inequality (14) holds, which is equivalent to 

-Xj(p) < Xi(p) < pXj(p). (20) 

l Ji 

Meanwhile, the closed-loop system satisfies the following LMI performance level 7 j over each pa- 
rameter subset Vi 

A d,i(p, p) x i(p) + Xi{p)AciMp) + J2{Pk^k} -rr- 1 + A 0 Xi(p) Xi{p)B cUl (p ) C^{p) 


BUp)Xi(p) 

C c l,i{p) 


- 7 -Jn d B d,i(p) 

B d,i{p ) — 7 i^n e 


where 


A cl,i{Pi P) B cl,i{p ) 

Ccl,i(p) D c pi{p) 


Mp) 

0 

BiAp) 

0 

B 2 Ap) 

0 

0 

0 + 

I 

0 

. c u (p) 

0 

DuAp) . 

0 

D 12 ,i(p) 


A k,i(PiP) B k,i(P) 0/0 

Ck,(p) D k ,i{p) C 2 ,i{p) 0 D 2 \,i{p) 

Then the closed-loop system is exponentially stable for every switching signal a with average dwell 
time (19) and its performance is maintained as ||e ||2 < 7 ||c /|| 2 with 7 = max { 7 i} ieZN - 

The proof of the exponential stability for average dwell time switching is similar to that in [7] and 
will be omitted here. In addition, given the initial condition x c ;(0) = 0, it can be shown from the 
condition ( 20 )— ( 21 ), the inequality 

K + -e T e - 7 d T d < 0 , 7 = max{~.;j ; 

7 

holds within each parameter subset. Integrate on both sides of the inequality, we get 


V a (x cl (T)) - V a ( Xcl (0)) + — ||e||| - 7 ||rf||l < 0. 

7 
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Since V a (x c i(T)) > 0 and V a (x c i( 0)) = 0, we have ||e ||2 < 7 ||d ||2 hold as desired. 

Using Elimination lemma, one can verify the equivalence between condition (21) and LMIs (15)- 
(17). Moreover, when Mi(p) = Ri(p) and Nj.(p) = R~ 1 (p) — Si(p), recall the expression of the 
Lyapunov function V a (x c i,p) in (9). If switching from the parameter subset V, to Vj, it is required 
that Vj([x;x\, p) < pVj(x c i. p) to guarantee the stability condition of switching with average dwell 
time. Obviously, it is satisfied if pRj(p) > Ri(p). Similarly, Ri{p) > jjRj(p) is sufficient to hold 
Vi([x;x\,p) < pVj(x c i,p) when switching from the other direction of the switching surface Sij. 
Combining them together gives the boundary condition (18). Q.E.D. 


Note that the switching LPV synthesis condition for this switching logic is different from hysteresis 
switching LPV control results. The (1,1) term in ( 15)— (16) implies that the closed-loop switched 
LPV system has its convergence rate at least The open- loop plant can be thought as a shifted 
system with its A,; matrix changing to A /. It is the same for controller A k ^ matrix. Therefore, 
if the matrix functions Ri{p) and Sj(p) can be solved, then the gains of the switching LPV controllers 
will be constructed by replacing A, and A k j in the standard LPV controller formula (10) by A, + 4/7 
and A k: , + \l. 

3 Flight Control Example 

The system to be controlled is the longitudinal F-16 aircraft model based on NASA Langley Re- 
search Center (LaRC) wind tunnel tests [24], which is described by Stevens and Lewis in great 
detail [25]. In this research, a simple thrust vectoring model is added to the aircraft model and 
will be activated in the high angle of attack region to provide additional longitudinal axis control 
power. 

3.1 Longitudinal Model of F-16 Aircraft with Thrust Vectoring 

The states used to describe the motion of the aircraft in longitudinal axis over entire operating 
envelope are as follows: V (ft/s) is the total aircraft velocity, a (deg) is the angle of attack, q 
(deg/s) is the pitch rate, and 6 (deg) is the pitch angle. For the original aircraft model, the 
available control inputs are the throttle setting 5th and the elevator angle 5 e (deg). The resulting 
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nonlinear equations of motion in longitudinal axis are given as follows: 


V = — (F x cos a + F* sin a) 
m 

(22) 

a = — — (— F x sin« + F „ cos a) + q 
mV 

(23) 

. My 

q =Y 

1 y 

(24) 

6 =q, 

(25) 


where m is the aircraft mass, F x and F z are the force components along x and z body axes 
respectively, I y is the moment of inertia about the y body axis, and M y is the pitching moment. 
Note that the throttle setting indirectly affects the states through the power output from the engine. 
Therefore, the actual power level is also considered as a state variable in longitudinal dynamics, 
and the detailed dynamical model of the engine can be referred in the NASA data. In addition, V. 
a, and q are selected as outputs. 

The x and z axes forces and pitching moment in Eqs. (22)-(25) contain aerodynamic, gravitational 
and thrust components. 

F x =qSC x j. — mg sin# + T x 
F z =qSC z j + mg cos 6 + T z 
M y =qScC m j + Mt, 

where q is the dynamic pressure, S is the wing surface area, and c is the wing mean aerodynamic 
chord. A complete description of the total coefficients C x x, C z j, and C m j can be founded in Ref. 
[24], which also provides the aerodynamic data in tabular form. 



(a) Original aircraft model 


(b) Aircraft model augmented with thrust vectoring 


Figure 3: Aircraft model with and without thrust vectoring. 


The F-16 aircraft is powered by an after-burning turbofan jet engine, which produces a thrust force 
in the x axis direction, as shown in Figure 3(a). In this research, the F-16 aircraft is augmented 
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with a simple thrust vectoring model, which is similar to that in Ref. [2], Denote the thrust vector 
angle by S p t v , as shown in Figure 3(b). With the right-handed (forward, starboard, and down) 
coordinate system, the thrust components along the x, z axes and the pitching moment due to 
thrust vector are then given by 


T x — T cos Spt v 
T z — T sindpf v 
Mt = — ItT sin< 5 p tv, 

where It is the moment arm from the center of gravity to the thrust application point. A more 
complicated model of thrust vectoring can be found in Refs. [26, 27], but will not be used in our 
study. 

To develop an LPV representation of the nonlinear F-16 model, we first need to find the wings- 
level equilibrium points at several flight conditions in the design envelope. The local linear models 
are then obtained by linearizing the nonlinear equations of motion at those points. The flight 
envelope of interest covers aircraft speeds between 160 ft/s and 200 ft/s and angles of attack from 
20° to 45°. These two variables are used as scheduling parameters in the LPV modeling of F-16 
longitudinal dynamics. The points at which the nonlinear model is linearized are marked by a “x” 
symbol in Figure 4. To apply the switching LPV control synthesis technique, the flight envelope is 
partitioned into two subregions. For hysteresis-based switching, the striped area in Figure 4(a) is 
the overlapped parameter region. So there are two switching surfaces, a = 30° and a = 35°. For 
the switching with average dwell time, the overlapped region shrinks to a line. Here, we take the 
line a = 33° shown in Figure 4(b) as the switching surface. 

In this research, two different sets of actuators are used in the different angle of attack regions. 
The actuators used in the low angle of attack region (region 1), are the throttle and the elevator, 
and the thrust vectoring nozzle is inactive. The local linear models in this region are based on 
the original aircraft model, which is corresponding to the case of S p t v = 0. In the high angle of 
attack region (region 2), the thrust vectoring nozzle is incorporated to provide additional force and 
moment. Therefore, the switching of controllers is based on the trajectory of the angle of attack, 
i.e., the controller is switched only when the aircraft flies from one angle of attack region to the 
other. 
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a (deg) a (deg) 

(a) Hysteresis switching (b) Switching with average dwell time 


Figure 4: Flight conditions and partitioned flight envelope for two switching logics. 

3.2 Control Problem Setup 

The control design objective is to track the command of the angle of attack. It is formulated as 
a model-following problem, where the ideal model to be followed is chosen to be a second-order 
low-pass filter based on desired flying qualities. A block diagram of the system interconnection 
for synthesizing the switched LPV controllers is shown in Figure 5, where P is the model set of 
linearized aircraft dynamics at different trim points, and it has three outputs: the velocity V , the 
angle of attack a, and the pitch rate q. The inputs of the open-loop system include 3-dimensional 
sensor noise signal n. the angle of attack command cn cm d, and the control input u. The outputs of 
the open-loop system are weighted error signals e p and e u , and the measurement y. 



Figure 5: Weighted open-loop interconnection for the F-16 aircraft. 

As mentioned before, pilots desire fast and accurate responses in the low angle of attack region. 
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While in the high angle of attack region, the requirement for the flying quality is not so critical, and 
the emphasis of the control design is the maintainability of the stability. Therefore, we choose the 
different performance weighting functions in the different angle of attack regions to reflect design 


requirements. 


TT , 100(5/20 + 1) w 80(5/10 + 1) 

W V\ ~ „ /n 1 1 1 ’ W p 2 ~ , nm , , 


P 5/0.1 + 1 ’ P 5/0.01 + 1 ’ 

where the subscripts 1 and 2 denotes the low and high angle of attack region, respectively. Note 
that the tracking error in the steady state is 1% in the region 1 and 1.25% in the region 2 . Also, 
the bandwidth in the region 1 is greater than that in region 2 , and this is consistent with the design 
objective. 


Also, different actuator sets are used in different angle of attack regions. The dynamics of the 
actuators are modeled as first order lag filters, and the time constants can be found in [24, 2], 

$th 1 dptv 1 $e 1 

“ 0.2.s + 1 ’ “ 0.07.S + 1 ’ “ 0.05.5 + 1 ' 

In the low angle of attack region, the control inputs are throttle position %/, and elevator angle S e . 

Both the positions and the rates of control inputs are fed into W u to penalize the control effort. 

Therefore, the system matrix of W act is given by 

-5 5 

W act = diag < 10, 

-5 5 

and the weighting function W u is given as 

W„ = diag{l,10,+ +}. 

In the high angle of attack region, the thrust vector is also activated. In order to use the proposed 
switching condition, the weighted open-loop LPV plants in the different parameter subspaces must 
have the same number of the states. Otherwise, the synthesis conditions on the switching surfaces 
will not hold. Therefore, we ignore the dynamics of the throttle, which has the slowest time constant 
among the three actuators, and keep the order of W act as same as that in the low angle of attack 
region. The related weighting functions W act and W u are given as follows. 

, W u -diagji, 35 , i 20 ’ 50 ’ 120 } ' 
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The other common weighting functions are chosen as 


W n = diag {0.8, 0.1, 0.6} , 


R^ideal 


6.25 

s 2 + 3.6s + 6.25’ 


where the ideal model is a second-order system with the natural frequency 2.5 rad/s, and the 
damping ratio 0.72. 


3.3 Design Results and Nonlinear Simulations 

Two switched LPV controllers corresponding to hysteresis-based and average dwell time switching 
logics are designed using the proposed switching LPV synthesis conditions in Theorem 1 and 2. The 
multiple parameter-dependent Lyapunov functions at each parameter subset are parameterized as 
affine functions of scheduling parameters. That is, we have for % = 1, 2 

Ri{p) = Ri + PiR} + P2R] t Si(p) = Si + PiS} + P2S ? , 

where p\ = a, p 2 = V, and matrices Rf and S-' : with A; = 0, 1, 2 are new optimization variables to 
be determined. The general objective function to optimize is defined as min^^ij wpp. where w, 
is the weight to penalize 7 ., and Ylf = 1 Ul i = 1- I 11 this paper, we define the objective function as 
minmaxjyi, 72 }, which means minimizing the “worst-case” performance level 7 ., . The performance 
level 7 over the whole parameter set are 7.9913 for hysteresis switching logic, and 15.9914 for 
average dwell time switching with Ao = 0.01 and p = 1 . 1 . 

For the switching with average dwell time, the influence of the constants Ao and p is also studied 
by comparing five different sets of values listed in Table 1. The data in the first row are those used 
in the previous control design and the next nonlinear simulation. We first keep Ao fixed at 0.01, 
and change the value of p. It is noticed that the variation of p just changes the value of the average 
dwell time t„. and has no impact on the performance level 7 . This can be explained by the fact 
that p is not involved in the synthesis conditions (15)— (18), and thus not related to the value of 7 . 
We then keep p = 1.1 and let Ao vary. The performance level 7 is improved if Ao decreases, but 
the resulted optimization problem for the case of Ao = 0.05 is infeasible. It is still an open question 
to determine suitable values of Ao and p. One possible way is to try and obtain an acceptable 
performance level as well as a reasonable average dwell time value. 

To test the performance of the switching control system during the nonlinear simulation, two angle 
of attack command inputs are defined as shown in Figure 6 . The initial angle of attack command is 
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Table 1: Effect of constants Aq and /i on average dwell time switching. 


Ao 

I 1 

r a 

7 

0.01 

1.1 

9.53 

15.9914 

0.01 

1.05 

4.88 

15.9914 

0.01 

1.5 

40.55 

15.9914 

0.001 

1.1 

95.31 

8.4199 

0.05 

1.1 

1.91 

infeasible 


selected as 36°, and thus the initial switching signal cr(0) = 2. The trajectory of the angle of attack 
command is defined as square waves, and deliberately chosen to cross the two parameter subsets 
back and forth to illustrate the effect of switching LPV control. There are four switching events 
happened. For the first command, the switches occur at 0 s, 15 s, 31 s, and 45 s, respectively. 
Therefore, all the time intervals between switches satisfy the requirement r a > 9.53 s. To see the 
effect of the average dwell time r a , we defined another command input which switches fast at the 
beginning and then compensates by switching sufficiently slow. The switches occur at 0 s, 5 s, 10 
s, and 45 s, respectively, and thus the average dwell time still satisfy the requirement. Note that 
the switching time of the actual angle of attack trajectory is slightly different from the command 
input, since we try to track the response of a second-order ideal model, not the square wave. 



40 
S 35 


30 

time (s) 


30 

time (s) 


40 


50 


(a) Command input 1 


(b) Command input 2 


Figure 6: Command inputs for switching LPV control simulation. 


The test flight condition is selected at V = 180 ft/s and a = 36°. Figure 7 shows the nonlinear 
responses of the aircraft model for the command input 1. The dotted lines in subplot (a) rep- 
resent the angle of attack response of the ideal model, the solid and dashed lines in all subplots 
represent the responses using hysteresis switching and average dwell time switching, respectively. 
The tracking performances over the entire time history are acceptable for both switching logics. 
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It is noticed that the switching signal a in subplot (b) is a little different for those two switching 
logics, since their switching surfaces are different. For hysteresis switching, the switch event occurs 
when a = 30° or 35°. While for average dwell time switching, there is only one switching surface 
a = 33°. The responses of the actuators are shown in subplots (d)-(f). Obviously, the thrust vector 
is activated only when the aircraft flies in the high angle of attack region. 



(a) Angle of attack 



(c) Velocity 



(e) Thrust vectoring nozzle deflection 
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Figure 7: Nonlinear simulation of switching LPY control for command input 1. 


Figure 8 presents the nonlinear response to the command input 2. For the second switching logic, 
although the first two switch events occur too fast, the average dwell time still satisfies the re- 
quirement (19). Therefore, the tracking performance of switching with average dwell time is also 
acceptable and a little bit worse than that achieved by the hysteresis switching. However, the 
stability of the switched system will not be guaranteed if the condition of the average dwell time 
is violated. Compared to the average dwell time switching, there is no such a restriction on the 
hysteresis switching. In addition, according to the theoretical developments and the simulation re- 
sults, the switching with average dwell time appears more complicated than the hysteresis switching 
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since one has to choose the parameter Ao and /t. Thus the hysteresis switching scheme seems more 
flexible and has a potential for the high performance aircraft applications. 



0 10 20 30 40 50 60 0 10 20 30 40 50 60 


time (s) time (s) 

(e) Thrust vectoring nozzle deflection (f) Elevator angle 

Figure 8: Nonlinear simulation of switching LPV control for command input 2. 

4 Conclusion 

A switching LPV control approach based on multiple parameter-dependent Lyapunov functions is 
presented for flight control design. The state of the controller is reset to the value of the plant 
state when the switching event occurs. The controller state reset not only guarantees the stability 
requirement of Lyapunov function, but also leads to the formulation of switching LPV control 
synthesis conditions as LMI optimization problems. Two parameter-dependent switching logics, 
hysteresis-based switching and switching with average dwell time, are used to avoid the possible 
transient instability caused by switching among controllers. The proposed switching LPV control 
technique is applied to an F-16 aircraft model switching between low and high angle of attack 


21 








regions with different control objectives and actuator sets, and promising simulation results are 
obtained. 
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